BCS-BEC Crossover in Bilayers of Cold Fermionic Polar Molecules 
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We investigate the quantum and thermal phase diagram of fermionic polar molecules loaded in a 
bilayer trapping potential with perpendicular dipole moment. We use both a BCS theory approach 
that is most realiable at weak-coupling and a strong-coupling approach that considers the two-body 
bound dimer states with one molecules in each layer as the relevant degree of freedom. The system 
ground state is a Bose-Einstein condensate (BEC) of dimer bound states in the low density limit 
and a paired superfluid (BCS) state in the high density limit. At zero temperature, the intralayer 
repulsion is found to broaden the regime of BCS-BEC crossover, and can potentially induce system 
collapse through the softening of roton excitations. The BCS theory and the strongly-coupled dimer 
picture yield similar predictions for the parameters of the crossover regime. The BKT transition 
temperature of the dimer superfluid is also calculated. The crossover can be driven by many-body 
effects and is strongly affected by the intralayer interaction which was ignored in previous studies. 

PACS numbers: 03.75.Ss,67.85.-d,74.78.-w 



I. INTRODUCTION 



Recent progress on trapping and cooling of polar 
molecules [lHSl enable studies of many-body systems 
with long-range anisotropic dipole-dipole forces, where 
new exotic phases could exist 0, Q ■ The attractive part 
of the interaction can, however, lead to chemical reac- 
tion losses [6]. A way to stabilize the system is to load 
molecules in a one- or two-dimensional optical lattices, 
where interesting low-dimensional physics has been pre- 
dicted 0-[T3|- A prominent example is polar fermions 
loaded in a bilayer system with dipoles oriented perpen- 
dicular to the layer plane [l8[ . In the limit of high density 
or weak interaction, the system is very similar to conven- 
tional superconductors, and the ground state should be 
a BCS state (with interactions as in Fig. [Ha)). In the 
dilute limit, it is known that the interlayer interaction 
always supports a bound state [191 - 121) . and the ground 
state should be a BEC of dimers (Fig.QJb)). As a result, 
the BCS-BEC crossover in this bilayer system could be 
richer than the usual atomic Fermi gas where crossover is 
driven by the two-body physics of a Feshbach resonance 
[221 ] . Here it is driven by many-body effects which depend 
not only on interaction strength but also on density. Fur- 
thermore, the intralayer repulsion can cause roton soft- 
ening and/or Wigner crystallization in the high density 
limit [15l4l7| |. Therefore, a many-body theory including 
both intra- and interlayer interactions is not only of quan- 
titatively interest but also qualitatively important. 

In this paper, we study the quantum and thermal 
phase diagrams of fermionic polar molecules loaded in 
a bilayer system including both intra- and inter-layer 



interaction. The quantum phase diagram is shown in 
Fig- Hie). BCS (BEC) ground states are found in 
the limit of weak (strong) interaction and large (small) 
density. In between, we have the crossover regime 
(Crossover-I), which can be determined by the chemi- 
cal potential calculated from many-body theories in the 
different limits (see below). When neglecting intralayer 
repulsion, the crossover region (Crossover-II) moves to 
lower interaction strength. In addition, intralayer in- 
teractions could perhaps give rise to a roton instability 
at large density and strong interaction, although further 
analysis beyond the scope of this work is needed. We 
also determine the BKT critical temperature [U [24| in 
the strongly-coupled regime including the effective inter- 
action of the dimers. The maximum critical temperature 
obtained is one tenth of the Fermi energy, and should 
therefore be achievable in experimental setups in the near 
future. 

We note that two other recent studies have considered 
a system similar to the one studied here. The paper by 
Pikovski et al. |25[ considers the BCS and BEC phases 



based on BCS theory at both zero and finite tempera- 
ture. However, they do not consider the full effect of the 
intralayer interaction. When we neglect the intralayer 
interaction, our results are consistent with Ref. [25) . A 
related study by Baranov et al. (26[ addresses the crit- 
ical temperature for the superfluid phase in the weak- 
coupling limit, taking also particle-hole correlations into 
account. In contrast, here we consider the finite tem- 
perature phase diagram from the strong-coupling limit. 
Extrapolation of our results to the parameter regime of 
Ref. |26| would exceed the boundaries of our approxima- 
tions and the current study should be viewed as comple- 
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FIG. 1: (color online) Schematic of a few particles in the 
bilayer with (a) inter- (Vi) and intralayer (Vb) interactions in 
the BCS limit and (b) the effective interaction (V e s) between 
the dimers in the BEC limit (see text), (c) Quantum phase 
diagram for a bilayer with fermionic polar molecules in the 
interaction-density plane. Interactions are characterized by 
U = mD 2 /h 2 d and by kpd, where m, D, d, and &f are re- 
spectively the molecular mass, dipole moment, interlayer dis- 
tance, and Fermi momentum. The Crossover-I region includes 
the effect of intralayer interactions, while the Crossover-II re- 
gion does not. In the upper right-hand corner we speculate 
that the system could potentially display a roton instability 
as discussed in the text. 



is the area of the layer plane; e^, = h 2 k 2 /2m, and 
a = ± is the layer index. V+_ = V-+ = V\ de- 
notes the interlayer interaction and V++ = V = Vo 

is the mirelayer one. Here we neglect interlayer tun- 
neling. We assume occupation of only the ground state 
in the transverse direction so that the transverse de- 
gree of freedom is a simple gaussian that can be inte- 
grated out. This yields V (q) = (i -\F{\q\)) 

and Vi(q) -> -2-KD 2 \q\e-^ d as W/d -> 0. Here q 
is the in-plane momentum, D is the dipole moment, 
W and d the layer width and interlayer spacing, and 
F(q) = s/n/2Wq[l - E^(Wq/V2)}e q2w ^ 2 with Erf(ir) 
the error function. The interlayer interaction is exact in 
the strict 2D limit (W <C d), but we find it is accurate 
enough (< 10%) at W = 0.2d, which is the value used 
throughout. In the weak interaction limit, the intralayer 
repulsion should renormalize the single particle disper- 
sion as in Fermi liquid theory. Treating the layer index 
as a spin degree of freedom, the gap equation is 
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where Eq = y£q + |Aq| 2 is the quasi-particle disper- 
sion, and (,q = eq + £(g) - hbcs- The self-energy, E(q), 
within in the Hartree-Fock approximation is [27| S(fe) 
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Eq(Vo(0)-V (k-q))[l-^tanh(E q /2k B T)]. To ; 



cess the crossover regime, the chemical potential, [ibcs, 
must be determined self-consistently via the density 
equation 
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mentary to Ref. [261 ] . 

Our model and the assumptions we use are described 
Section [TTJ In Section Hill we discuss the case of zero tem- 
perature and the role of intralayer interactions on the 
BCS-BEC crossover. This is achieved by considering the 
physics from both a weak- and a strong-coupling point 
of view. Both approaches are shown to yield consistent 
results. We proceed to discuss the finite temperature 
phase diagram in Section HVl Here we calculate the crit- 
ical temperature for the superfluid phase within mean- 
field theory and using the universal relation for the BKT 
transition temperature. Section [V] contains a summary, 
a discussion of experimental parameters to realize the 
predicted phases, and an outlook for future work. 



II. MODEL 

The Hamiltonian for polar molecules in the bi- 
layer system is given by H = J2ka e fc c fc. CT c fc, CT + 

£iEqkk\aa> V ™'^ c k+qJk' -qy C k' ,o> c k^ where n 



We use the first Born approximation (FBA) with a real- 
istic finite layer width for both intra- and interlayer inter- 
actions in Eq. ([T]). The exponentially decreasing shape 
of V\{q) means that using the renormalized gap equa- 
tions or the FBA yields similar results as noted already 
in Ref. [25(. The FBA is generally poor at strong cou- 
pling. However, since the crossover takes place at low 
density, we treat the intralayer in a weak coupling sense 
in the current work. Inclusion of higher-order terms will 
be considered in future work. 

Since the interlayer interaction, V\, is attrative at 
short distance, the dominant gap is s-wave. Therefore, 
for simplicity, we neglect higher partial wave compo- 
nents and calculate A^, = A,^, from Eq. ([T]) by iter- 
ation. In Fig. [H[a), we show Afe/Ao for different val- 
ues of the dimensionless coupling, U = mD 2 /h 2 d with 
kpd = V Annd 2 = 0.4. When U is small, the maximum 
of A(fc) is at k/kp — 2.5, which is where V\(k) is most 
attractive. When U is larger, the intralayer interaction 
renormalizes the single particle dispersion; the structure 
of A/,, disappears and the maximum is at k — 0. In 
Fig.[i;b) we show A (T)/E F as function of T/T F , where 
Ep and Tp are Fermi energy and temperature. 
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U=l k F d = 0.4 




FIG. 2: (color online) (a) Normalized gap, Afc/Ao, as a func- 
tion of fc/fc_F at kpd = 0.4 and {/ = 1 (solid) and [7 = 5 
(dashed), (b) Temperature dependence of Ao(T)/ Ef as a 
function of T/Tf for fc F d = 0.4 = 0.16£ ). (c) BCS 

critical temperature as a function of kpd at U = 2 in units of 
So = h 2 /2md 2 . (d) Ratio of the roton wavelength, X r , at the 
instability to the dimer size, Ib, as function of kpd. 



III. ZERO TEMPERATURE RESULTS: 
INTRALAYER INTERACTION EFFECTS ON 
THE BCS-BEC CROSSOVER 

When the density is reduced, the ground state is ex- 
pected to be a dimcr BEC. Within BCS theory, this 
regime can be defined by having a negative chemical po- 
tential {pLscs < 0). From Eq. (0), we can determine 
the hbcs — boundary as shown in Fig. [He). The two 
solid lines that bound the two crossover regions are cal- 
culated by including (fiBes) and neglecting (^bcs) ^ ne 
intralayer repulsion respectively. As expected, the repul- 
sion strongly suppresses the hbcs < region and for 
densities larger than a critical value of kpd ~ 0.77, /j-bcs 
is positive for any U. We always find a ^%cs < ® region 
when the intralayer interaction is neglected. This shows 
that the intralayer interaction brings not only quanti- 
tative contributions, but also qualitative and important 
changes of the quantum phase diagram. The effect of 
intralayer repulsion on the crossover physics can also 
be seen in the transition temperature, T c , as shown in 
Fig. Hfc). It has a maximum at fcp-d ~ 1, similar to 
crossover in a quasi-2D superconductor (28[. In contrast 
to BCS results for Fermi gases with short-range interac- 
tion, the decrease of T c at high density is caused by the 
long-range intralayer repulsion. 

In order to investigate the many-body physics of the 
BEC limit (i.e. the strong interaction or low density 
limit) in more detail, we go beyond BCS theory, which 
is based on Fermi liquid theory at high density. Starting 
from the extremely dilute limit, the bound dimer is the 
main constituent with binding energy Eb- For small den- 
sity, we neglect the intralayer term but include the Fermi 



pressure in the chemical potential; uPp ^ n = Ep — Eb/2, 
where Eb is obtained numerically 19 -[21 [ . We define the 
crossover regime (Crossover-II) as the region bounded by 
Hbcs = and ^bec = 0; the latter given by the line 
with triangles in Fig. [He). The Crossover-II region is 
quite narrow which indicates that the BCS theory and 
the strong-coupling result are very similar. The expres- 
sion for Mnr- g g iven above holds exactly for zero-range 
interactions [29[, and indicates that the crossover hap- 
pens when the dimer size becomes comparable to the 
interparticlc distance. While the interlayer dipole inter- 
action behaves similar to a zero-range interaction, the 
intralayer dimer-dimer interaction can have significant 
effects on (Ibec as we now demonstrate by deriving an 
effective interaction between dimers. We note that a re- 
cent Monte Carlo study of the BCS-BEC crossover in two 
dimensions also finds that in the BEC limit, the dimer- 
dimer and atom-dimer interactions are important correc- 
tions that are not taken into account in the usual BCS 
theory without self-energy corrections [13] ■ 

The coordinates of the four molecules are denoted by 
ri, • • • ,T4 as shown in Fig. [TJb), where (r±,r2) are for 
the left dimer and {r^.r^) for the right dimer. We are 
interested in deep bound state, where the dimer size is 
smaller than the inter-dimer distance, i.e. \p\ » \r\, |r'|, 
where p = (ri + t*2)/2 — (r^ + r*4)/2 is inter-dimer dis- 
tance, and r = r\—T2 and r' = r% — are relative coor- 
dinates in each dimer. Straighforward algebra and inte- 
gration over the dimer bound state wavefunction, <p(r), 
yields the effective dimer-dimer interaction V e g(p) = 
J drdr' |0(r)| 2 |0(r')| 2 E s =±i, Q =o,i V a {p + sr a ) 7 where 
ro,i = (r =p r')/2. In the strong interaction and di- 
lute density limit, we can approximate 4>{r) by a Gaus- 
sian profile; <p(r) = {IbV^)^ 1 cxp (— |r| 2 /2Z 2 3 ), where 
l B = ^/2h 2 /niE B is the radius of dimcr bound state. As 
a result, V e g has the following simple Fourier transform: 

V cS (k) = [2Vx(k) + 2V (k)] exp(-\k\ 2 l 2 B /8). (3) 

The Gaussian approximation for </>(r) fails for U < 2, 
however we have checked that our results are qualita- 
tively unchanged if the exact solution is used [3l| . No- 
tice that V e g takes the strong interlayer interaction into 
account through the dimer wavefunction. The FBA is 
used for the intralayer term, which is reasonable since 
the crossover happens at low density. For the roton in- 
stability, the FBA was used in quasi-2D studies of dipolar 
bosons [32T [33| and we expect it to remain a fair approx- 
imation at larger densities as well. We note that dipolar 
interactions are different from short- or zero-range in- 
teractions. Zero-range interactions in 2D always allow 
a two-body dimer bound state [29[, and a state of four 
bosons will also be bound [Hj]. This result is different 
from 3D where the interaction must be sufficiently at- 
tractive to produce bound states. In the current setup, 
the dimer-dimer system is unbound [3oT |. 

The effective interaction has the property that V c s(0) 
is non-zero due to the intralayer term. However, the 
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molecules are fermionic and the Fock exchange contri- 
bution could cancel this term as is the case with true 
short-range interactions in interacting single-component 
Fermi gases. However, V e g(q) also has large contribu- 
tions from qd > 0, and this should have an influence on 
the phase diagram. This is supported by recent studies of 
the density- wave instability where the effect of exchange 
is found to be very large, shifting the instability into the 
strong-coupling regime p3a - l40l |. The importance of ex- 
change effects has also been discussed in relation to fer- 
roeletricity with polar molecules [4lj| . Here we estimate 
the effects in the bosonic dimer limit by including the 
intralayer potential through V c g(0). 

Using V c g , the chemical potential in the BEC limit can 
be estimated as [lbec = ^Kff(0)/2 — E B /2 (we neglect 
the Fermi pressure which is much smaller than the in- 
teraction energy for strong interactions). In Fig. Hfc), 
the solid line with crosses is given by /j-bec — and is 
the upper bound of the crossover regime including in- 
tralayer interaction (Crossover-I). Notice in particular 
that the BEC region shrinks to lower density when in- 
cluding the intralayer interaction, and there is no dimer 
condensate for kpd > 0.42. The ratio of interparticle dis- 
tance to dimer size is five or more along the (J-bcs = and 
/J-bec = lines. Without intralayer interaction the ra- 
tio is around 2.5, again demonstrating how the crossover 
physics is strongly modified by the intralayer term. We 
have estimate the chemical potential from both BCS and 
BEC limits, giving bounds of the crossover regime with 
and without the intralayer interaction. Our results in- 
dicate that in a realistic experiment, the intralayer re- 
pulsion can significantly affect the regime where a dimer 
condensate is observable. 

From the effective interaction between dimers, V e g, 
we can also calculate the dispersion, Tuujd(k) — 
^/e^,/2(e^,/2 + 2nV c g (A;)), of the collective Bogoliubov 
mode of a dimer supcrfluid in analogy to the case of dipo- 
lar bosons [33]. With increasing U we find roton soften- 
ing around kd = 2tt, which corresponds to a wavelength 
of X r = 2-K/k ~ d. This softening leads to system col- 
lapse in the high density and strong interaction regime, 
as shown in Fig. [TJc). To investigate the nature of the 
instability, we plot \ t /Ib in Fig. HJd) and observe that 
A r is more than a factor of two larger than Is for all den- 
sities. This implies that, at least in the low-density limit, 
it is a many-body effect. 

The roton analysis assumes a well-defined dimer pic- 
ture. We find the instability in the upper right-hand 
corner of Fig. [Tt, i.e. at higher densities and large U. 
While the dimer picture should prevail for large U, the 
higher densities imply that some fermionic nature could 
perhaps arise. Similar effects should arise in a fermionic 
picture but it is not easily calculated since fluctuations 
beyond BCS theory are needed. So at this point the find- 
ings for the roton instability remains speculative. This is 
an interesting topic for future work. We note that den- 
sity waves in a single layer with fermionic polar molecules 
have been predicted recently [ID, |43[ . A rough estimate 



indicates that the roton instability lies inside the region 
where a density wave in single layers is predicted. How- 
ever, Refs. J42j and [43[ neglect the Fock contributions 
which are expected to be large [36M40I ] . The exact region 
of the density wave instability is therefore not yet known. 
In any case, we expect the system to be unstable in the 
upper right part of the zero temperature phase diagram 
of Fig. □ 



IV. FINITE TEMPERATURE RESULTS: 
INTRALAYER INTERACTION EFFECTS ON 
THE CRITICAL TEMPERATURE 

We now investigate the finite temperature phase di- 
agram. In the BCS limit, we can use Eq. ((T|) to ob- 
tain the transition temperature, T Cl as shown by the 
dashed line in Fig. [3] We note that the intralayer in- 
teraction has very little influence on T c at kpd — 0.4. It 
is known that in the weak interaction limit, T c is very 
close to the true transition temperature in a full BKT 
theory [Hj]. However, this calculation fails in the strong 
interaction or dilute limit. At strong coupling, we ob- 
tain the BKT transition temperature from the univeral 
relation ksTsKT — fi 2 "Kn s (TBKT) /2m, where n s (T) is 
the superfluid density at temperature T [H, H3- Ac- 
cording to Landau's two-fluid model, we have n s (T) = 
n — n n (T), with normal fluid density given by n n (T) — 

h 2 f d 2 q 
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In Fig. [3] we 



[Sinh(uj d (q)/2k B T)_ 

show Tbkt of a dimer superfluid as a function of U for 
kpd = 0.4 (solid black line). The dimer result differs 
from the BCS theory prediction in the weak coupling 
limit, while it becomes saturated at Tbkt = 0.125Tp 
for large U. In Ref. (25|, T c and Tbkt was calculated 
from the BCS superfluid state. For small U we find the 
same result for T c . For dipolar bosons, Ref. [46[ finds 
Tbkt ~ O.lTp for all values of U shown in Fig. [3] In 
order to compare to that study we need to assume point 
dimers with twice the dipole moment, which is of course 
not the case for small U where Tbkt is reduced as the 
dimer size grows. Our strong-coupling dimer approach 
to the finite temperature physics can be considered com- 
plementary to both Ref. [25[ and Ref. (46|. The good 
agreement at intermediate and large U with the other ap- 
proaches provides support that we capture the essential 
physics using the effective interaction between dimers. 



V. DISCUSSION 

Wc have studied fermionic polar molecules in a bilayer 
with perpendicularly polarized dipole moments. As the 
density and dipole strength varies, an analogue (but with 
different physics) of the celebrated BCS-BEC crossover 
is predicted. We find that the intralayer repulsion (and 
dimer-dimer interaction) is crucial quantitatively and 
also to some extend qualitatively: It shrinks the region 
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FIG. 3: (color online) Transition temperature (Tbkt, solid 
black line) as a function of interaction strength for fc^d = 0.4. 
BCS result (T c , dashed blue line) is shown for comparison. 
The dashed vertical line marks the position of (ibcs = (see 

Fig. mc)). 



of dimcr condensation so that no dimer state is expected 
when the density is larger than a critical value and it also 
causes a roton instability in the high density and strong 
interaction regime. Our work is thus important for the 
study of BCS-BEC crossover physics in 2D systems, and 
our results should be observable within the parameter 
regime of near future experiments. 

It is worth stressing that our study considers the 
physics of the bilayer with fermionic polar molecules from 
different points of view. The BCS theory is usally more 
reliable in the weak-coupling regime, but can be extended 
into the crossover regime by solving the gap and number 
equations self-consistently. As the bilayer setup will al- 
ways have a two-body bound state with one molecule in 
each layer, it is reasonable to consider these dimers as 
the relevant degree of freedom in the strongly-coupled 
regime. The results presented above indicate that both 
approaches yield consistent results, with or without in- 
cluding the repulsive intralayer interaction. 

The crystalline phases [15l - [l7| are also ground state 
candidates due to the intralayer repulsion. For W/d = 
0.2 used in our calculations, these phases could appear 
below the roton instability region in Fig[TJ The finite ex- 
tend of the dimers is, however, still a concern and further 



work is needed to determine crystal phases in the bilayer 
setup. 

In order to detect the phases a number of techniques 
could be applied. Dimerization in the layers should 
be detectable by Bragg scattering [47| or in situ non- 
demolition detection [48|, [49| , whereas RF spectroscopy 
can probe the gap. To probe the finite temperature 
physics one can detect the associated vortices by matter 
wave heteredyning [50] . To estimate parameters for rele- 
vant systems, we take d = 0.5/im which yields U = 1.22 
for KRb molecules at D = 0.566 Debye and U = 4.15 for 
LiCs at D = 1.0 Debye (the permanent dipole moment 
of this molecule in the ground-state is about D = 5.4 
Debye [51j). For densities n = 10 5 — 10 8 cm -2 we 
have kpd = 0.06 — 1.8. The critical density to reach 
Mscs < is about n ~ 0.2 • 10 8 cm" 2 , whereas the 
dimer BEC requires n < 0.6 • 10 7 cm" 2 . Unfortunately, 
as Tp = E^kpd) 2 /ks this means that extremely low 
temperatures (< InK) are required to reach Tbkt- 

Interesting directions for future work includes more 
than two layers or tilting of the dipoles with respect to 
the plane. With dipoles that are no longer perpendicu- 
lar one can still show that a two-body bound dimer will 
be present for any value of the dipole moment [5^ „ al- 
though the dimer binding energy is reduced [53| . In the 
single layer case, p-wave superfluidity can occur for an 
extended range of tilting angles [ll| and we expect simi- 
lar effects for a bilayer. This setup can be explored with 
the same methods used here. In a setting with multiple 
layers more phases should be expected such as coherent 
density waves (36|, pairin g (471 ■ and bound states with 
more than two molecules 13 5l I54J1 . 
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